Finite-Size Scaling in the Energy-Entropy Plane for the 2D ±J Ising Spin Glass 
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For L x L square lattices with L < 20 the 2D Ising spin glass with +1 and -1 bonds is found 
to have a strong correlation between the energy and the entropy of its ground states. A fit to the 
data gives the result that each additional broken bond in the ground state of a particular sample of 
random bonds increases the ground state degeneracy by approximately a factor of 10/3. For x = 0.5 
(where x is the fraction of negative bonds), over this range of L, the characteristic entropy defined 
by the energy-entropy correlation scales with size as L 1 ' 78 ^ 2 ' . Anomalous scaling is not found for 
the characteristic energy, which essentially scales as L 2 . When x — 0.25, a crossover to L 2 scaling 
of the entropy is seen near L — 12. The results found here suggest a natural mechanism for the 
unusual behavior of the low temperature specific heat of this model, and illustrate the dangers of 
extrapolating from small L. 



I. INTRODUCTION 

The Edwards- Anderson (EA) spin glassi has been 
studied extensively for thirty years. A complete under- 
standing of its behavior in two and three dimensions re- 
mains elusive. In recent years it has become possible 
to compute the free energy of the two-dimensional (2D) 
Ising spin glass with ± J bonds on L x L lattices with L 
of 100 or morei^iMiSiSii From these calculations on large 
lattices we have learned that extrapolations of data from 
lattices with L < 30 are often misleading i&2ii2*ii*i2i 

A better understanding of why this happens is clearly 
desirable. This is especially true because essentially all 
of the work on three-dimensional (3D) EA models at low 
temperatures must be done on lattices with L < 20, due 
to our inability to equilibrate larger lattices at low tem- 
peratures in 3D>i£ At least one example of complex be- 
havior of the order parameter emerging as L is increased 
is already known in a similar 3D model ^* 

In this work we will analyze data for the energies and 
entropies of the ground states (GS) of 2D Ising spin 
glasses obtained using methods from earlier worki&2ii£ 
We will demonstrate that for small square lattices the 
± J EA model has a strong correlation of the sample-to- 
sample fluctuations of the energy and the entropy of the 
GS. The increase of GS entropy So with GS energy E is 
too large to be explained by fluctuations in the number 
of zero-energy single-spin flips. This correlation may be 
the cause of the breakdown of naive scaling behavior at 
small L in this model. 



II. THE MODEL 



The Hamiltonian of the EA model for Ising spins is 



H — ^ JijSiSj , 

m 



(i) 



where each spin Si is a dynamical variable which has 
two allowed states, +1 and -1. The (ij) indicates a sum 



over nearest neighbors on a simple square lattice of size 
L x L. We choose each bond Jy to be an independent 
identically distributed quenched random variable, with 
the probability distribution 

P(Jij) = xS{J i3 + 1) + (1 - x)5(J ij ~ 1) , (2) 

so that we actually set J = 1, as usual. 

The data analyzed here used an ensemble in which, for 
a given value of x, every L x L random lattice sample 
had exactly (1 — x)L 2 positive bonds and xL 2 negative 
bonds. Details of the methods used to calculate Eq and 
the numbers of GS have been described earlier^ 5 Sq is 
defined as the natural logarithm of the number of ground 
states. For each sample, Eq and Sq were calculated for 
the four combinations of periodic (P) and antiperiodic 
(A) toroidal boundary conditions along each of the two 
axes of the square latticed We will refer to these as 
PP, PA, AP and AA. We use ALL to refer to a data set 
which includes the results from all four types of boundary 
conditions. In the spin-glass region of the phase diagram, 
the variation of the sample properties for changes of the 
boundary conditions is small compared to the variation 
between different samples of the same sizejS except when 
x is close to the ferromagnetic phase boundary and the 
ferromagnetic correlation length becomes comparable to 
L. 



III. GROUND STATE PROPERTIES 

The average GS entropy (So) of an L x L sample for 
this model is essentially proportional to L 2 , the number 
of spins, with a small finite-size correction^ It was dis- 
covered earlier^ however, that for x = 0.5 the ratio of 
the width of the distribution of Sq for different samples 
of size L divided by (Sq) is not a monotonic function of 
L, having a peak at L — 8. A similar change in behavior 
between L — 8 and L = 10 was seen earlier by Saul and 
Kardar in samples with open boundary conditions, and 
appears in Fig. 11 of their paper & The original motiva- 
tion of the current study was to understand the origin of 
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FIG. 1: (color online) Scatter plots of correlations for x = 0.5 and L = 10: (a) Eo averaged over boundary conditions vs. 
number of frustrated plaquettes; (b) So averaged over boundary conditions vs. number of frustrated plaquettes; (c) So vs. Eo; 
(d) Eq vs. So- The number of samples used is 400, and the lines through the data are least-squares fits. 



this unexpected behavior. Our ensemble, unlike the one 
used by Saul and Kardar, does not keep the number of 
frustrated plaquettes fixed. 

We first look to see if the GS properties are correlated 
with the number of frustrated plaquettes, with the num- 
ber of bonds of each type held fixed. The scatter-plot 
data for x — 0.5 and L = 10 arc shown in Fig. 1(a). 
There is a substantial correlation of Eq with the number 
of frustrated plaquettes, and this correlation seems to be 
independent of L. Since it is well known that Eq increases 
as the number of frustrated plaquettes is increased, this 
is expected. 

There is a weaker correlation between So and the num- 
ber of frustrated plaquettes, as shown in Fig. 1(b). On 
the average, increasing the number of frustrated plaque- 
ttes increases So. This correlation is also not surprising, 
since positive So arises from rearranging the strings of 
broken bonds which connect the frustrated plaquettes in 
a GS. It seems natural that a larger number of frustrated 
plaquettes would give a larger number of ways to rear- 



range the strings of broken bonds. 

For Fig. 1(a) and 1(b), we have averaged Eq and Sq 
over the four different boundary conditions for each sam- 
ple, because the number of frustrated plaquettes does not 
depend on the boundary conditions. In the remainder of 
this work, we will treat each boundary condition for each 
sample independently. 

All equilibrium statistical mechanics can be derived 
from the partition function, which is determined by the 
energy and the entropy. Therefore, we would like to know 
if Eq and Sq are correlated with each other. The scatter 
plots for this correlation from the same data are shown 
in Fig. 1(c), along with a least-squares fit to the data, 
treating Eq as the independent variable. Fig. 1(d) shows 
the same data with the roles of energy and entropy re- 
versed. It demonstrates that the least-squares fit depends 
on which variable is chosen as the independent one. 

The results of least-squares fits for x — 0.5, 0.25 and 
0.125, and L varying from 6 to 20 are shown in Figures 
2, 3 and 4, respectively. For each value of x and each L, 
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FIG. 2: (color online) Results of least-squares fit analysis pa- 
rameterized by Eqn. (3) of the scatter plot of correlations 
between Eo and So for x = 0.5: (a)b vs. L, log-log plot; 
inset: b/L 2 vs. L; (b)m vs. L. The number of samples 
used for each L, (L, #), is (6:400), (8:400), (10:400), (12:400), 
(14:400), (16:400), (18:400) and (20:238). 



we show the slope m given by the least-squares fit, and 
the offset b of the entropy, defined by 

S'(E ) = m*E + b. (3) 

S' is the dependent variable in the least-squares fit. Note 
that Eo is negative. 

In Figures 5, 6 and 7, we give the value of r, the nor- 
malized covariance defined by 

r{E ,S ) = — , (4 

a{Eo)a(So) 

for each fit. The angle brackets indicate an average over 
the random bond distribution for some fixed value of x. 
The standard deviation, a, is defined, as usual, as 



a(X) = VV 2 )- (X) 2 . (5) 
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FIG. 3: (color online) Results of least-squares fit analysis pa- 
rameterized by Eqn. (3) of the scatter plot of correlations 
between Eo and So for x = 0.25: (a)6 vs. L, log-log plot; 
inset: b/L 2 vs. L; (b)m vs. L. The number of samples 
used for each L, (L, #), is (6:200), (8:200), (10:200), (12:200), 
(14:200), (16:200), (18:133) and (20:200). 



The numbers in parentheses and the error bars shown in 
the figures represent a one standard deviation statisti- 
cal error, as calculated by the Origin 6.0 Professional 17 
least-squares fitting routine. One expects that, in addi- 
tion, there may be systematic errors arising from non- 
ideal behavior of random number generators and nonlin- 
ear correlations. It is often difficult to obtain meaningful 
estimates of systematic errors. 

For small L there is a strong correlation between Eq 
and So. As L increases, a(Eo) increases linearly with L 
but <j(So) increases faster than linearly over this range 
of L. Thus the correlation gets weaker as L increases. 
This is reflected in the tendency for r to decrease as L 
increases. From our data, it is not clear whether or not r 
goes to zero as L goes to infinity. It is generally believed 
that the model is not self-averaging at T — 0, so it would 
be natural for r to remain finite as L increases. 

The value of r depends on the choice of ensemble. For 
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FIG. 4: (color online) Results of least-squares fit analysis pa- 0.4 -| — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , — , 
rameterized by Eqn. (3) of the scatter plot of correlations 
between Eo and So for x — 0.125: (a)6 vs. L, log-log plot; 
inset: b/L 2 vs. L\ (b)m vs. L. The number of samples 
used for each L, (L,#), is (6:200), (8:200), (10:200), (12:200), 

(14:200), (16:200), (18:200) and (20:200). FIG. 7: (color online) Covariance of ground state energy and 

ground state entropy, r(Eo,So) vs. L for x = 0.125. 
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FIG. 5: (color online) Covariance of ground state energy and 
ground state entropy, r(Eo, So) vs. L for x = 0.5. 



the ensemble in which we fix both the number of nega- 
tive bonds and the number of frustrated plaquettes for 
each value of L, one would find higher values of r than 
what we find here. Crudely, one would expect the values 
of r in the more tightly specified ensemble to be higher 
by about a factor of 1/.8, the inverse of the r-factor for 
the correlation between Eq and the number of frustrated 
plaquettes. 

For x = 0.25 and x = 0.5, as L increases the slope of 
the regression line through the data given by the least- 
squares fit appears to rapidly approach a limit of about 
TTi w 0.36. This number is slightly greater than Zn(2)/2 = 
0.34657.... Naively, this means is that, on the average, 
the GS degeneracy increases by about a factor of two 
for each additional broken bond, since each broken bond 
increases the energy by two units. But m is not actually 
a physical observable, because it depends on our choice 
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FIG. 8: (color online) Results of least-squares fit analysis pa- 
rameterized by Eqn. (6) of the scatter plot of correlations 
between Eo and So for x — 0.5: (a)— A vs. L, log-log plot; 
inset: -A/L 2 vs. L; (b)B vs. L. 



FIG. 9: (color online) Results of least-squares fit analysis pa- 
rameterized by Eqn. (6) of the scatter plot of correlations 
between Eo and So for x = 0.25: (a)— A vs. L, log-log plot; 
inset: -A/L 2 vs. L; (b)B vs. L. 



of ensemble. We will say more about what the actual 
physical quantity is later. 

The reader should note that the probability den- 
sity in the energy-entropy plane shown in Fig. 1(c) is 
clearly different from a two-dimensional Gaussian distri- 
bution, since, with the boundary conditions we are us- 
ing, Eq can only have values which are multiples of four 
units. The previous results in Figure 6 of Landry and 
Coppersmith^ using a much larger number of samples, 
show that the one-dimensional probability distribution 
for So at the same values of x and L, which is the projec- 
tion of the joint distribution onto the entropy axis, can 
be fit by a Gaussian distribution. 

Since it is generally believed that T = is a critical 
point for the model, perhaps what one needs to explain is 
why the one-dimensional distribution is apparently Gaus- 
sian! This result becomes less surprising, however, once 
one realizes that the strong correlations between Eq and 
the number of frustrated plaquettes will make it very 



difficult to see any non-Gaussian behavior in this one- 
dimensional distribution, unless one holds the number of 
frustrated plaquettes fixed. In our ensemble, the number 
of frustrated plaquettes always has a Gaussian distribu- 
tion. 

For x — 0.125, where the ferromagnetic correlations 
are substantial for small L, the strength of the energy- 
entropy correlation is somewhat reduced for the case of 
periodic boundary conditions in both directions. This 
effect is probably a result of the fact that for small L 
at x = 0.125, the behavior is essentially dominated by 
short-range ferromagnetic correlations. 

It is equally valid to do the least-squares fit using So as 
the independent variable. The results of fits of this type 
for x = 0.5, 0.25 and 0.125, using the same data as before, 
are shown in Figures 8, 9 and 10, respectively. The values 
of r are not shown again, since they are unchanged from 
the earlier case. For each value of x and each L, we now 
show the slope B of the least-squares fit, and the offset 
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FIG. 10: (color online) Results of least-squares fit analysis 
parameterized by Eqn. (6) of the scatter plot of correlations 
between Eo and So for x = 0.125: (a)— A vs. L, log-log plot; 
inset: -A/L 2 vs. L; (b)B vs. L. 



A of the energy, defined by 

E'(S ) = A + B*S 



(6) 



Now the notation E' indicates that energy is the depen- 
dent variable in the least-squares fit. Since A is negative, 
we use —A in the log-log plot. 

From the log- log plots in Figs. 2(a) and 3(a), we find 
that b is proportional to L 178 ( 2 ) and L 181 ' 4 ', respec- 
tively, while from Figs. 8(a) and 9(a), the scaling of —A 
is essentially indistinguishable from L 2 . A closer inspec- 
tion of Fig. 3(a), however, reveals that there is a clear 
curvature of the data on the log-log plot. We can see 
in the inset to Fig. 3(a) that, for x = 0.25, b becomes 
essentially proportional to L 2 for L > 12. 

There is no known reason for a qualitative difference 
in the scaling behavior between x = 0.5 and x = 0.25. 
Therefore, we anticipate that the entropy scaling will also 
become proportional to L? for some larger L in the x — 
0.5 case. Given the results in the literature, we expect 
that this will happen before L reaches 30. 

As we have demonstrated in Fig. 1(c) and 1(d), the 
regression line which is obtained when one uses So as 
the independent variable is not the same one which is 
found by using Eq as the independent variable. The slope 
parameters of these two regression lines are related to 
each other— as 



mB 



(7) 



Since r 2 is, in general, a number between and 1, we 
are faced with the problem of deciding what the true 
best line through the data is. Without some additional 
information, there is no unique prescription for solving 
this problem^ 

We can write down an equation for the joint probabil- 
ity distribution which builds in the fact that the allowed 
values of Eq are quantized: 



P L {E ,S )=C L f( 



E - (E (L)) S Q -(S (L)) . 
a(E (L)) ' a(So(L)) ■ 



E S(Eo-4n)] 



(8) 



r 



The dependence on x is not shown explicitly, and Cl is 
the normalization constant. Aside from small corrections 
to scaling which can be ignored for large L, we expect 
that we can assume 



(E (L)) = E^L 2 



(9) 



and 



(S (L)) = SooL 2 . (10) 
It should also be sa foV 5 to assume that ct(Eq(L)) = 



oeL- The scaling of <j(Sq(L)) with L appears to be 
nontrivialj2ii£ but we certainly expect that it will diverge 
as L goes to infinity. Therefore, for large L it should be 
an adequate approximation to replace the sum over the 
(5-functions by a uniform background. 



If the envelope function / of the probability distribu- 
tion in the energy-entropy plane was a two-dimensional 
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Gaussian, then it would have the normal form 



fc(X,Y) 



27TV1 



: exp 



(- 



X 2 - 2rXY + Y 2 
2(1 -r 2 ) 



(11) 

where the arguments X and Y have probability distribu- 
tions with zero mean and unit standard deviation. Given 
this assumption, which is unproven, we should treat Eq 
and So on an equal basis. Then it would be correct to 
set the best regression line through the joint probability 
distribution to be equal to 



So — S 
or, equivalcntly, 



(m/r) * (Eb - £oc) , 



E -E oo = (B/r) * (S - Soo) . 



(12) 



(13) 



As we have repeatedly reminded the reader, T = is 
believed to be a critical point for this model, so the as- 
sumption of Gaussian fluctuations can be justified only 
as an approximation. We do not really know what ac- 
tual form of f(X, Y) should be used, and therefor e) 19 ! 20 
we do not know what the slope of the best regression line 
should be. However, it seems certain that corrections to 
the Gaussian approximation are invisible at our current 
level of statistical uncertainty. 



IV. DISCUSSION 

One might think that the behavior of the A and B pa- 
rameters which we find by treating the entropy as the in- 
dependent variable in the least-squares fit are quite sim- 
ple. It is important to remember, however, that the slope 
of the best line through the data is not B. Within the 
Gaussian approximation, as we have remarked above, the 
slope of the best line is B/r. And thus if the slope of the 
best line through the data is to have some finite slope 
in the large L limit, it appears necessary to have a finite 
limit for r. 

In Figures 11, 12 and 13, we show the values of B/r 
vs. L. The value of B/r appears to be approximately 
independent of L for x = 0.125 and 0.25, because the L 
dependences of B and r cancel each other, although the 
statistical uncertainties are too large for precise state- 
ments to be made. For x = 0.5, however, B/r may be 
increasing with L within our range of L. 

Dimensionally, the slope B/r defined in Eqn. (13) has 
units of temperature. It is tempting to argue that B/r 
has some relation to a fictive glass temperature for the 
crossover between high and low temperature dynamical 
behavior. Thus, a naive prediction for B/r would be that 
it should be proportional to the mean-field energy scale, 
which is 



E, 



mf 



2y/x(l-x), 



(14) 



where the factor of 2 comes from the square root of the 
number of neighbors on the lattice, and the factor of 
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FIG. 11: B/r vs. L for x = 0.5. 
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FIG. 12: B/r vs. L for x = 0.25. 



x(l — x) is the second moment of P(Jy). However, no 
such dependence on x is seen in our data. The value of 
B /r actually seems to be decreasing slowly as x increases 
from 0.125 to 0.5. For x = 0.5, if we average the data for 
all L, we find 



B/r = 1.66 ±0.03. 



(15) 



This number probably underestimates the result for large 
L slightly, due to the apparent tendency for B/r to in- 
crease with L. The quoted statistical error does not in- 
clude any allowance for this effect. Using Eqn. (15), how- 
ever, we find that, on the average, each additional broken 
bond in the ground state increases the GS degeneracy by 
a factor of 3.34(7). Our uncertainties for the smaller val- 
ues of x are larger, but this is primarily because we have 
smaller numbers of samples for these cases. 

For any GS on a square lattice, each spin which has 
two neighbor spins that are oriented along the direction 
favored by the bond between them, with the other two 
neighbor spins pointed opposite to the direction favored 
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FIG. 13: B/r vs. L for x = 0.125. 



by the bond (i.e. these bonds are broken), can flip with 
no energy cost. Each of these free spins contributes a 
factor of two to the degeneracy of the GS. It is thus ex- 
pected that increasing the number of broken bonds in the 
GS would also increase the number of such free spins. It 
is not reasonable, however, that increasing the number 
of broken bonds by one would increase the number of 
free spins by nearly two, on the average. Therefore, the 
fact that we have found the average increase in the GS 
degeneracy for each additional broken bond to be a fac- 
tor of about 10/3 indicates that this effect cannot be 
explained by fluctuations in the number of zero-energy 
single-spin flips. There must be a substantial contribu- 
tion from large-scale rearrangements of the GS structure. 

An exact calculation of the energy-entropy correlation 
for L = 50, or possibly L = 60, could be performed using 
the method of Galluccio, Loebl and VondrakiSi It should 
be noted, however, that it is not really necessary to calcu- 
late So exactly. It would be more than sufficient to have 
an approximate calculation of Sq which was accurate to 
one part in 10 4 . That does not seem impossible, and it 
might allow calculations for even larger L. 

An explicit calculation of the low temperature specific 
heat for x — 0.5 by Lukic et al& gives a result propor- 
tional to exp(— 2/T) when L > 30. It is natural that 
crossover behavior should be seen in both 5*0 and the low 
T specific heat, with the same crossover length. The low 
T limit of the specific heat in this model is determined 
by the degeneracy of the states at energy 4 above the 
GS. It is surely not surprising that the degeneracies of 
the ground states and the first excited states would be 
controlled by the same crossover length. To verify that 
this is occurring, the calculation of Lukic et al. could be 
repeated for x = 0.25. We expect that a crossover length 
of L w 12 will be found for the specific heat in that case. 

Katzgraber and Lee 2 ^ have calculated the T depen- 
dence of the correlation length in this model, and found 
that it behaves as exp(2/T). (Recall that J = 1.) They 



use this result to argue that the specific heat at low 
T should be proportional to exp(— 4/T), as one might 
naively expect for a model with an energy gap of 4. How- 
ever, a more detailed analysis^ has found that their data 
for the specific heat agree with the conclusions of Lukic 
et al& Another recent study by WangjSi using a new 
algorithm, also finds that the low T specific heat is pro- 
portional to exp(— 2/T). 

It would also be interesting to repeat these calcula- 
tions on a hexagonal lattice, where the allowed energy 
states are multiples of two units, because the number of 
bonds for each site is odd. Thus on this lattice the small- 
est zero-energy excitation of the ± J model involves two 
neighboring spins. We would expect that the low tem- 
perature specific heat is proportional to cxp(— 2/T) for 
the whole range of L in that case. 

The analogy to an Ising chain which is made by Wang 
and Swendseni^ to argue for a specific heat which is pro- 
portional to exp(— 2/T) has nothing to do with random 
bonds. We know, however, that in 2D a fully frustrated 
Ising system does not display this behavior. 9 In addition, 
by studying triangular lattices, Poulter and Blackman 25 
found that adding a small concentration of unfrustrated 
plaquettes to a fully frustrated system does not produce 
spin-glass behavior. 

It was recently shown by Amoruso, Marinari, Martin 
and Pagnani 26 that the behavior of domain wall energies 
for the 2D Ising spin glass is fundamentally different in 
those cases, such as the ± J model, where the energies are 
quantized. In a very interesting paper, Wang, Harrington 
and PreskillSi have argued that the presence of an energy 
gap allows the existence of topological long-range order. 
A less specific suggestion of topological long-range order 
in 2D random-bond Ising models was made earlier by 
Merz and Chalkcr. 28 Numerical results for the properties 
of domain walls, which will be presented elsewhere (SLSl 
are consistent with this proposal. 

What we want to do is to explain the low temperature 
behavior of the specific heat in terms of topological ex- 
citations. We know that in the spin- glass region of the 
phase diagram it is difficult to find a way to overturn a 
finite fraction of a large sample at zero energy cost. In 
contrast, at and near the fully frustrated system, where 
the exp(— 2/T) behavior of the specific heat does not oc- 
cur, it becomes easy to find ways of overturning a finite 
fraction of the spins at zero energy cost. This explains 
the difference between the spin-glass region and the fully 
frustrated region. 

What is possible in the spin-glass region is finding 
many ways of overturning a finite fraction of the large 
sample at a cost of only 4i2S This is precisely the general- 
ization of the Wang-Swendsen mechanism 16 to a random 
2D system. And because we can do this, the behavior of 
the Parisi overlap function- a will be nontrivial. Of course, 
since T c is zero, the overlap function collapses to zero as 
L increases. But if we scale out this simple collapse with 
L, the existence of the large-scale finite-energy excita- 
tions may be observed in the overlap function. Since we 
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are trying to observe effects caused by states of energy 
4 J, we cannot merely do a naive T — calculation of the 
overlap function. However, by manipulating the double 
limit L — > oo and T — ► 0, it may be possible to see the 
effect. 

Although the actual implementation would be very 
challenging, one can imagine studying the energy-entropy 
correlation at and below T c in a three-dimensional Ising 
spin glass, using thermal-average values for the energy 
and entropy. In that case, where the spin-glass phase and 
the failure of self-averaging of the Parisi overlap function 
are believed to occur at finite one can use a gen- 
eral probability distribution for P(Jij) and still have a 
positive entropy, in contrast to the 2D situation. 

The author's expectation is that the anomalous scaling 
which we find for the ± J model at T = in 2D for small 
L will occur for all types of bond distributions in 3D, 
where the spin-glass transition is at T > 0. It might also 
happen that the crossover length becomes infinite in 3D, 
but it seems more likely that the crossover length is only 
infinite in four or more dimensions. 



V. SUMMARY 

We have found that for L x L square lattices with 
L < 20 the 2D Ising spin glass with +1 and -1 bonds 
has a very strong correlation between Eq and the number 
of frustrated plaquettes, and, what is more surprising, a 
strong correlation between Eq and So- On the average, 



each additional broken bond in the GS of a particular 
sample of random bonds increases the GS degeneracy 
by a factor of about 10/3. This number is too large to 
be explained by fluctuations in the number of free spins, 
which implies that there is a substantial contribution due 
to large-scale rearrangements of the GS structure. Over 
this range of L, the characteristic GS entropy scales as 
£!-78(2) £ or x _ o.5 5 while the characteristic GS energy 
scales cis , <xs expected. For x = 0.25, however, a 
crossover is seen to normal scaling behavior of Sq near 
L = 12. We believe that a similar crossover will occur for 
x = 0.5 at L » 25, and that this crossover is connected 
to the anomalous behavior of the low temperature spe- 
cific heat. We explain why the Wang-Swendsen mecha- 
nism for a low T specific heat proportional to exp(— 2/T) 
should apply in the spin-glass regime, but not in the fully 
frustrated regime. 
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